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In August 201 1, an anthrax outbreak occurred among Hippopotamus amphibius hippopotamuses and humans in Zambia. Here, 
we report the draft genome sequence of the Bacillus anthracis outbreak strain CZC5, isolated from tissues of H. amphibius hip- 
popotamuses that had died in the outbreak area. 
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Anthrax is a life-threatening disease caused by the Gram- 
positive spore-forming bacterium Bacillus anthracis (1). The 
disease occurs worldwide, and outbreaks and epidemics of the 
disease in both humans and animals have been reported in Zam- 
bia (2, 3). In August 2011, an anthrax outbreak occurred among 
Hippopotamus amphibius hippopotamuses in the Chama district 
of the eastern province of Zambia (4). Over 80 hippopotamuses 
died after showing signs of infection with B. anthracis. Following 
the deaths of the hippopotamuses, 521 suspected human cases 
resulting in 6 deaths were reported. As of January 2014, 7 com- 
pleted genome sequences and 23 draft genome sequences for 
B. anthracis had been deposited in GenBank. However, most of 
these belong to strains that were isolated in the United States; there 
have been no reports on the genome sequences of B. anthracis 
strains isolated in central Africa. Here, we report the draft genome 
sequence of a B. anthracis strain associated with an outbreak 
linked to unnatural deaths of hippopotamuses in Zambia in 201 1. 

The genomic DNA was extracted from overnight cultures us- 
ing a Qiagen DNeasy blood and tissue kit, according to the sup- 
plied protocol. A library for sequencing was prepared from the 
extracted genomic DNA using a Nextera XT sample preparation 
kit (Illumina), and paired-end sequencing with a read length of 
2 X 300 bp was performed on an Illumina MiSeq platform. The 
total amount of read data comprises 1,223,422,928 bases, repre- 
senting 234-fold coverage of the genome. After obtaining reads 
that passed the quality threshold with the Trimmomatic-0.20 tool 
(5) and merging overlapping paired-end reads with Fastqjoin 
(https://code.google.eom/p/ea-utils/), a mapping assembly was 
carried out with Newbler version 2.9 (Roche). The assembled ge- 
nome consists of 25 contigs of >301 bp (N 50 contig size, 
195,537 bp), with an average G+C content of 35.12%. The draft 
genome was structurally and functionally annotated using the 
Rapid Annotations using Subsystems Technology (RAST) anno- 
tation server (6). A total of 5,841 putative coding sequences 
(CDS), 1 1 rRNA operons, and 91 tRNA genes were identified. The 



quality-passed reads were aligned to the reference B. anthracis 
strain Ames Ancestor (taxonomy ID 261594) using SSAHA2 (7) as 
part of the Breseq pipeline version 0.24rc4 (8), resulting in only 1 
missing gene (that for aspartate ammonia-lyase [aspA-3]) in this 
strain. 

We believe that a detailed comparative genomic analysis of this 
B. anthracis strain, designated CZC5, with more genome se- 
quences of this species obtained in Zambia might lead to new 
insights into its evolutionary history and pathogenicity. The re- 
sults of such studies combined with data for worldwide lineages 
will be included in a future publication. 

Nucleotide sequence accession number. This whole-genome 
shotgun project for B. anthracis strain CZC5 has been deposited at 
DDBJ/ENA/GenBank under the accession no. BAVT00000000. 
The version described in this paper is the first version. 
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